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ISSUE:  Dredged  material  disposal  sites  located 
in  open  water  are  classified  as  either  dispersive 
or  nondispersive  depending  on  whether  local 
water  velocities  are  strong  enough  to  erode  and 
transport  dredged  material  from  the  deposited 
mound.  The  Corps  needs  the  capability  to  predict 
stability  of  the  mound  and  long-term  migration 
patterns  of  eroded  material  to  (1)  identify  accept¬ 
able  disposal-site  locations,  and  (2)  provide  a 
quantitative  approach  for  gaining  site-designation 
approval. 

RESEARCH:  The  overall  work-unit  objective 
is  development  of  a  systematic  approach  for 
predicting  the  dispersion  characteristics  of  a 
specific  open-water  disposal  site.  This  objec¬ 
tive  includes  the  following  goals: 

•  Identify  realistic  wind-,  wave-,  tide-,  and 
storm-generated  velocity  boundary  condi¬ 
tions. 

•  Develop  numerical  models  capable  of  sim¬ 
ulating  dispersion  characteristics  of 
dredged- material  mounds  for  periods  of 
time  in  excess  of  one  year. 

•  Provide  site-designation  technology  to 
field  engineers  as  a  tool  in  site  identifica¬ 
tion  and  designation. 


This  study  reports  the  tidal-constituent 
database  developed  for  locations  along  the 
East  and  Gulf  coast  of  the  United  States, 
which  is  a  part  of  the  first  two  goals  listed 
above. 

SUMMARY:  The  numerical  model  ADCIRC- 
2DDI  (a  two-dimensional  depth-integrated  finite- 
element- based  hydrodynamic  circulation  code) 
was  applied  to  the  western  North  Atlantic,  Gulf 
of  Mexico,  and  Canibbean  Sea  to  develop  a 
tidal  constituent  database.  The  report  presents 
initial  results  from  a  tidal  model  of  the  western 
North  Atlantic  (WNAT)  that  encompasses  the 
coastal  ocean  as  well  as  the  deep  ocean.  The 
domain  for  the  model  includes  the  entire  area  of 
interest. 
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Preface 


The  research  described  in  this  report  was  authorized  and  funded  under 
Work  Unit  32466,  “Numerical  Simulation  Techniques  for  Evaluating  Long- 
Term  Fate  and  Stability  of  Dredged  Material  Disposed  in  Open  Water,”  of 
Technical  Area  1  (TA1),  “Analysis  of  Dredged  Material  Placed  in  Open 
Water,”  of  the  Dredging  Research  Program  (DRP).  The  DRP  is  sponsored 
by  Headquarters,  U.S.  Army  Corps  of  Engineers  (HQUSACE).  Adminis¬ 
trative  responsibility  is  assigned  to  the  U.S.  Army  Engineer  Waterways  Ex¬ 
periment  Station  (WES)  Coastal  Engineering  Research  Center  (CERC), 

Dr.  James  R.  Houston,  Director,  and  Mr.  Charles  C.  Calhoun,  Jr.,  Assistant 
Director. 

Mr.  John  H.  Lockhart,  Jr.,  HQUSACE,  was  DRPTA1  Technical  Moni¬ 
tor.  Mr.  E.  Clark  McNair,  Jr.,  was  DRP  Program  Manager  (PM),  and 
Dr.  Lyndell  Z.  Hales  was  assistant  PM.  Dr.  Nicholas  C.  Kraus,  Senior  Re¬ 
search  Scientist,  CERC,  was  the  TA1  Technical  Area  Manager,  and 
Dr.  Norman  W.  Scheffner,  Coastal  Processes  Branch  (CPB),  Research  Di¬ 
vision  (RD),  CERC,  was  the  Principal  Investigator  for  the  work  unit. 

Dr.  Scheffner  worked  under  the  administrative  supervision  of  Mr.  H.  Lee 
Butler,  Chief,  RD,  CERC,  and  Mr.  Bruce  A.  Ebersole,  Chief,  CPB,  RD, 
CERC. 

This  study  was  performed  and  the  report  was  prepared  over  the  period 
1  September  1990  through  15  December  1992.  The  numerical  modeling 
goals,  concepts,  and  methodology  were  developed  by  Dr.  Scheffner  and 
two  DRP  contractors:  Dr.  Joannes  J.  Westerink  of  the  University  of  Notre 
Dame  and  Dr.  Richard  A.  Luettich,  Jr.,  of  the  University  of  North  Carolina 
at  Chapel  Hill.  Development  and  implementation  of  the  model  were  com¬ 
pleted  by  Drs.  Westerink  and  Luettich. 

At  the  time  of  publication  of  this  report.  Director  of  WES  was  Dr.  Rob¬ 
ert  W.  Whalin.  Commander  was  COL  Leonard  G.  Hassell,  EN. 

Additional  information  on  this  report  can  be  obtained  from 
Mr.  E.  Clark  McNair,  Jr.,  DRP  Program  Manager,  at  (601)  634-2070  or 
Dr.  Norman  W.  Scheffner,  Principal  Investigator,  at  (601)  634-3220. 


VI 


Summary 


This  report  describes  the  application  of  model  ADCIRC-2DDI,  a  two- 
dimensional,  depth-integrated,  finite-element-based  hydrodynamic  circula¬ 
tion  code,  to  the  western  North  Atlantic,  Gulf  of  Mexico,  and  Caribbean 
Sea  in  order  to  develop  a  tidal  constituent  database.  Issues  that  are  empha¬ 
sized  in  the  development  of  the  Western  North  Atlantic  Tidal  (WNAT) 
model  include  definition  of  hydrodynamically  simple  open  ocean  bound¬ 
aries;  the  use  of  large  domains;  the  importance  of  a  high  degree  of  grid 
resolution  in  coastal  regions;  and  the  use  of  finite  element  meshes  with 
highly  varying  nodal  densities  in  order  to  minimize  the  size  of  the  discrete 
problem.  The  development  of  an  optimal  graded  finite-element  mesh  is 
based  on  regular  and  graded  grid  convergence  studies  using  an  M2  tidal 
forcing  function  on  the  boundary  and  within  the  domain. 

The  optimal  graded  mesh  is  then  forced  for  eight  diurnal  and  semidiurnal 
astronomical  tidal  constituents  0\,  P\,  Q\,  M2,  S2,  N2,  and  K2)  on  the 
open  ocean  boundary  by  coupling  to  Schwiderski ’s  (1979;  1981  a-g) 
global  model  results  as  well  as  within  the  interior  domain  using  a  tidal  po¬ 
tential  forcing  function.  Structures  of  the  various  tides  are  examined  and 
results  are  compared  to  field  data  at  77  stations. 
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1  Introduction 


There  has  been  a  recent  trend  in  coastal  ocean  tidal  modeling  towards 
using  increasingly  larger  computational  domains,  which  extend  up  to  or 
beyond  the  continental  shelf  break  and  slope.  For  example,  Flather  (1987), 
Gerritsen  and  Bijlsma  (1988),  and  Vincent  and  Le  Provost  (1988)  have  all 
developed  tidal  and/or  storm  surge  models  of  a  large  portion  of  the  eastern 
North  Atlantic.  These  studies  indicate  that  accurate  tidal  predictions  can 
be  conveniently  obtained  using  large  computational  domains.  This  report 
presents  initial  results  from  a  tidal  model  of  the  western  North  Atlantic, 
which  encompasses  the  coastal  ocean  as  well  as  the  deep  ocean.  Key  fea¬ 
tures  of  the  Western  North  Atlantic  Tidal  (WNAT)  model  are  the  defini¬ 
tion  of  hydrodynamically  simple  open  ocean  boundaries  to  facilitate  the 
specification  of  boundary  conditions,  the  use  of  a  high  degree  of  selective 
grid  refinement  to  resolve  the  flow  physics  on  a  localized  basis,  and  the 
coupling  of  the  coastal  ocean  with  a  global  tidal  model. 

The  domain  for  the  WNAT  model  encompasses  the  Western  North  At¬ 
lantic  Ocean,  the  Caribbean  Sea,  and  the  Gulf  of  Mexico.  The  WNAT 
model  has  an  eastern  open  ocean  boundary  that  lies  along  the  60°  W  merid¬ 
ian  and  is  situated  almost  entirely  in  the  deep  ocean.  The  location  of  this 
open  ocean  boundary  was  specifically  selected  to  simplify  the  difficult 
task  of  specifying  an  accurate  set  of  boundary  conditions  for  a  coastal 
ocean  tidal  model.  The  WNAT  model  open  ocean  boundary  offers  a  vari¬ 
ety  of  significant  advantages.  First,  the  boundary  is  geometrically  simple 
and  includes  no  discontinuities  or  corners.  Second,  the  boundary  lies  al¬ 
most  entirely  in  the  deep  ocean  where  tides  vary  gradually.  This  avoids 
the  difficulty  of  specifying  complex  and  highly  variable  cross  shelf  bound¬ 
ary  conditions  over  most  of  the  open  ocean  boundary.  Furthermore,  this 
open  ocean  boundary  is  ideally  suited  for  coupling  with  global  tidal  models, 
which  should  be  most  accurate  in  the  deep  ocean.  Third,  nonlinear  tidal 
constituents  will  not  be  significant  on  this  open  ocean  boundary  since  they 
are  generated  on  the  continental  shelf  and  are  largely  trapped  on  the  shelf 
due  to  the  out-of-phase  reflective  character  of  the  continental  slope. 

The  size  of  the  resulting  computational  domain  of  the  WNAT  model  ex¬ 
ceeds  8  x  106  km2  and  is  therefore  very  large  in  terms  of  coastal  ocean 
tidal  models.  It  is  clear  that  uniformly  discretizing  the  entire  domain  with 
the  resolution  required  in  regions  of  rapidly  varying  flow  would  substantially 


Chapter  1  Introduction 


2  Governing  Equations  and 
Numerical  Discretization 


The  computations  described  in  this  report  were  performed  using  ADCIRC- 
2DDI,  the  depth-integrated  option  of  a  system  of  two-  and  three-dimen¬ 
sional  hydrodynamic  codes  named  ADCIRC  (Luettich,  Westerink,  and 
Scheffncr  1992).  ADCIRC-2DDI  uses  the  depth-integrated  equations  of 
mass  and  momentum  conservation,  subject  to  the  incompressibility, 
Boussinesq,  and  hydrostatic  pressure  approximations.  Using  the  standard 
quadratic  parameterization  for  bottom  stress  and  neglecting  baroclinic 
terms  and  lateral  diffusion/dispersion  effects  leads  to  the  following  set  of 
conservation  statements  in  primitive  non-conservative  form  expressed  in  a 
spherical  coordinate  system  (Flather  1988,  Kolar  et  al.  1992): 
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t,  +  h  =  total  water  column 

bathymetric  depth  relative  to  the  geoid 

2Q  sin<p  =  Coriolis  parameter 

angular  speed  of  the  earth 

atmospheric  pressure  at  the  free  surface 

acceleration  due  to  gravity 

effective  Newtonian  equilibrium  tide  potential 

reference  density  of  water 

applied  free  surface  stress 
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Cj  =  bottom  friction  coefficient 


A  practical  expression  for  the  effective  Newtonian  equilibrium  tide  po¬ 
tential  is  given  by  Reid  (1990)  as: 


T)  (fc,  <p.  0  =  X  ajn  Cjn  fjn  W  Lj  (9)  C°S 


n,j 


2n  (t  - 10) 


Tjn+ft  +  vjn  ('(,). 


where 


Cjn  =  constant  characterizing  the  amplitude  of  tidal 
constituent  n  of  species  j 

otjn  =  effective  earth  elasticity  factor  for  tidal 
constituent  n  of  species  j 


fjn  =  time-dependent  nodal  factor 


v  •  =  time-dependent  astronomical  argument 


j  =  0,1,2  =  tidal  species  ( j  =  0,  declinational;  j  -  1,  diurnal; 
j  =  2,  semidiurnal) 


(4) 
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Lq  =  3  sin2  <p  -  1 
L,  =  sin  (2  <p) 

L2  =  cos2  (<p) 

X,  cp  =  degrees  longitude  and  latitude,  respectively 
t0  =  reference  time 

Tjn  =  period  of  constituent  n  of  species  j 

Values  for  Cjn  are  presented  by  Reid  (1990).  The  value  for  the  effective 
earth  elasticity  factor  is  typically  taken  as  0.69  for  all  tidal  constituents 
(Schwiderski  1980,  Hendershott  1981)  although  its  value  has  been  shown 
to  be  slightly  constituent  dependent  (Wahr  1981,  Woodworth  1990). 

To  facilitate  an  FE  solution  to  Equations  1-3,  these  equations  are 
mapped  from  spherical  form  into  a  rectilinear  coordinate  system  using  a 
Carte  Parallelogrammatique  (CP)  projection  (Pearson  1990): 


x'  =  R(X-X^)  cos  <p0 


(5) 


/  =  *<P 


(6) 


where  X.Q  ,  <p0  =  center  point  of  the  projection.  Applying  the  CP  projec¬ 
tion  to  Equations  1-3  gives  the  shallow-water  equations  in  primitive  non¬ 
conservative  form  expressed  in  the  CP  coordinate  system: 
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Utilizing  the  FE  method  to  resolve  the  spatial  dependence  in  the  shallow- 
water  equations  in  their  primitive  form  gives  inaccurate  solutions  with  se¬ 
vere  artificial  near  2  •  Ar  modes  (Gray  1982).  However,  reformulating 
the  primitive  equations  into  a  GWCE  form  gives  highly  accurate,  noise- 
free,  FE-based  solutions  to  the  shallow-water  equations  (Lynch  and  Gray 
1979,  Kinnmark  1984).  The  GWCE  is  derived  by  combining  a  time-differ¬ 
entiated  form  of  the  primitive  continuity  equation  and  a  spatially  differen¬ 
tiated  form  of  the  primitive  momentum  equations  recast  into  conservative 
form,  reformulating  the  convective  terms  into  non-conservative  form, 
and  adding  the  primitive  form  of  the  continuity  equation  multiplied  by  a 
constant  in  time  and  space,  Xq  (Lynch  and  Gray  1979;  Kinnmark  1984; 
Luettich,  Westerink,  and  Scheffner  1992).  The  GWCE  in  the  CP  coordi¬ 
nate  system  is: 
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The  GWCE  (Equation  10)  is  solved  in  conjunction  with  the  primitive  mo¬ 
mentum  equations  in  non-conservative  form  (Equations  8  and  9). 
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The  high  accuracy  of  GWCE-based  FE  solutions  is  a  result  of  their  ex¬ 
cellent  numerical  amplitude  and  phase  propagation  characteristics.  In 
fact,  Fourier  analysis  indicates  that  in  constant  depth  water  and  using  lin¬ 
ear  interpolation,  a  linear  tidal  wave  resolved  with  25  nodes  per  wave¬ 
length  is  more  than  adequately  resolved  over  the  range  of  Courant 
numbers  (C  =  'JghAt/Ax  <  1.0  (Luettich,  Westerink,  and  Scheffner 
1992)).  Furthermore,  the  monotonic  dispersion  behavior  of  GWCE-based 
FE  solutions  avoids  generating  artificial  near  2  •  Ax  modes,  which  plague 
primitive-based  FE  solutions  (Platzman  1981,  Foreman  1983).  The  mono¬ 
tonic  dispersion  behavior  of  GWCE-based  FE  solutions  is  very  similar  to 
that  associated  with  staggered  finite  difference  solutions  to  the  primitive 
shallow-water  equations  (Westerink  and  Gray  1991).  GWCE-based  FE  so¬ 
lutions  to  the  shallow-water  equations  allow  for  extremely  flexible  spatial 
discretizations,  which  result  in  a  highly  effective  minimization  of  the  dis¬ 
crete  size  of  any  problem  (Le  Provost  and  Vincent  1986,  Foreman  1988, 
Vincent  and  Le  Provost  1988,  Westerink  et  al.  1992). 

The  details  of  ADCIRC,  the  implementation  of  the  GWCE-based  solu¬ 
tion  to  the  shallow-water  equations,  are  described  by  Luettich,  Westerink, 
and  Scheffner  (1992).  As  most  GWCE-based  FE  codes,  ADCIRC  applies 
three-noded  linear  triangles  for  surface  elevation,  velocity,  and  depth.  Fur¬ 
thermore,  the  decoupling  of  the  time  and  space  discrete  form  of  the  GWCE 
and  momentum  equations,  time-independent  and/or  tri-diagonal  system 
matrices,  elimination  of  spatial  integration  procedures  during  time-stepping, 
and  full  vectorization  of  all  major  loops  results  in  a  highly  efficient  code. 
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3  Description  of  the 
Computational  Domain 


The  WNAT  model  domain  is  shown  in  Figure  1.  The  open  ocean 
boundary  extends  from  the  vicinity  of  Glace  Bay  in  Nova  Scotia,  Canada 
to  the  vicinity  of  Corocora  Island  in  eastern  Venezuela.  The  computational 
domain  is  bounded  on  the  north,  west,  and  south  by  the  North,  Central, 
and  South  American  coastlines.  The  topography  within  the  domain, 
shown  in  Figure  2,  includes  the  continental  shelf  (depths  typically  range 
from  0  m  to  130  m),  the  continental  slope  (depths  typically  range  from 
130  m  to  3,000  m),  and  the  continental  rise  and  deep  ocean  (depths  from 
approximately  3,000  m  to  almost  8,500  m). 

The  WNAT  model  open  ocean  boundary  lies  almost  entirely  in  the  deep 
ocean,  well  away  from  the  continental  shelf  and  slope.  As  was  discussed 
in  the  introduction,  it  is  highly  advantageous  to  specify  boundary  condi¬ 
tions  in  deep  water  since  tides  vary  slowly,  global  tidal  models  tend  to  be 
more  accurate,  and  nonlinear  constituents  are  small  in  deep  water.  In 
water  where  the  model  open  ocean  boundary  does  cross  or  come  near  the 
continental  shelf  or  slope,  there  are  nearby  islands  with  long-term  tidal 
data  which  allow  for  adjustments  if  the  global  model  data  prove  to  be  in¬ 
accurate  in  these  shallow  waters.  For  example,  the  very  northern  portion 
of  the  model  open  ocean  boundary  crosses  the  continental  shelf  and  slope 
but  also  intersects  Sable  Island,  Canada,  which  has  available  tidal  data. 
Furthermore,  the  southern  portion  of  the  open  ocean  boundary  nears  or 
crosses  the  continental  shelf  and  slope  but  also  is  in  the  vicinity  of  the 
lesser  Antilles  Islands  for  which  long-term  constituent  data  exist. 

Within  the  WNAT  model  domain,  147  stations  have  been  defined  at 
which  tidal  elevation  time  histories  were  recorded  for  all  model  runs.  These 
stations,  shown  in  Figure  2,  provide  a  basis  for  inter-grid  comparisons 
and/or  comparisons  with  field  measured  data. 

Stations  1  through  77  have  field-measured  long-term  tidal  elevation  data 
available  in  harmonically  decomposed  constituent  form  (long-term  indi¬ 
cates  that  the  time  history  records  used  to  determine  the  constituent  data 
are  typically  a  year  in  duration).  Details  regarding  the  locations  of  the 
stations  and  source  of  the  tidal  data  are  listed  in  Table  1 .  Data  at  a  given 
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station  from  the  various  sources  generally  agree.  These  tidal  elevation  sta¬ 
tions  lie  predominantly  in  open  waters  near  the  coast  along  the  eastern 
U.S.  seaboard  (stations  1-20),  within  the  Gulf  of  Mexico  (stations  21-38), 
and  within  the  Caribbean  Sea  (stations  39-58).  The  remaining  stations 
(stations  59-77)  are  generally  located  away  from  tne  coast  either  near 
small  islands  or  in  deep  waters  throughout  the  WNAT  model  domain.  Sta¬ 
tions  1  through  77  were  used  in  the  grid  convergence  studies  as  well  as  in 
comparing  the  computed  tides  with  the  measured  field  data. 

In  addition  to  stations  with  field-measured  tidal  data,  an  additional  70 
stations  were  defined  throughout  the  domain  in  order  to  achieve  a  more  ho¬ 
mogeneous  distribution  of  stations  within  the  domain.  The  additional  sta¬ 
tions  were  placed  on  the  continental  shelf  (stations  78-87),  in  the  vicinity 
of  the  continental  shelf  break  (stations  88-109),  near  the  toe  of  the  conti¬ 
nental  slope  (stations  110-130),  and  in  the  deep  ocean  (stations  131-147). 
These  stations  provide  a  more  thorough  basis  for  intercomparisons  of  re¬ 
sponses  for  the  grid  convergence  studies. 
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4  Grid  Convergence  Studies 


In  this  section,  the  resolution  requirements  for  the  WNAT-model  for 
the  semidiurnal  astronomical  M2  tide  are  examined.  A  sequence  of  three 
numerical  experiments  was  conducted  to  establish  the  resolution  required 
to  obtain  a  converged  M2  response  on  a  regular  grid,  to  develop  an  opti¬ 
mal  graded  grid  which  minimized  the  number  of  nodes  within  the  domain 
yet  gave  an  M2  solution  that  was  equivalent  to  the  finest  regular  grid  and 
to  examine  the  effect  of  resolution  of  the  boundary  on  the  M2  solution. 

The  three  numerical  experiments  involved  five  regular  grids  (R 1  -  R5)  and 
three  graded  grids  (VI  -  V3).  The  grid  properties  are  summarized  in 
Table  2.  The  coarsest  and  finest  regular  grids  are  shown  in  Figures  3  and 
4,  while  the  three  graded  meshes  are  shown  in  Figures  5  through  7.  All 
grids  were  generated  using  GREDIT,  a  flexible  interactive  grid  generation 
code  developed  by  Turner  and  Baptista  (1991). 

The  bathymetry  of  each  grid  was  obtained  by  interpolating  values  from 
the  ETOP05  database  from  the  National  Center  for  Atmospheric  Research, 
which  has  a  resolution  of  5'  by  5'.  However,  the  actual  representation  of 
the  bathymetry  for  each  grid  depends  on  the  grid  itself,  since  depth  is  in¬ 
terpolated  linearly  on  the  three-noded  triangular  finite  elements.  A  speci¬ 
fied  minimum  bathymetry  between  4  m  (for  the  Caribbean  Sea  and  Gulf  of 
Mexico)  and  7  m  (for  the  northeast  U.S.  and  Canadian  coasts)  was  speci¬ 
fied  to  avoid  the  drying  of  elements. 

The  open  ocean  boundary  for  each  grid  was  forced  with  an  M2  tidal  am¬ 
plitude  and  phase  interpolated  onto  the  open  ocean  boundary  modes  using 
data  from  Schwiderski’s  (1979)  global  model  results.  This  results  in  slightly 
different  boundary  forcings  in  space  for  grids  Rl,  R2,  and  R3.  Grids  R4, 
R5,  and  VI  through  V3  have  essentially  identical  boundary  forcing  func¬ 
tions,  An  effective  tidal  potential  forcing  was  applied  within  the  domain 
for  the  M2  constituent  using  an  amplitude  of  0.242334  m  (Reid  1990).  A 
theoretical  value  for  the  effective  earth  elasticity  factor  (which  reduces 
the  magnitude  of  the  tidal  potential  forcing  due  to  earth  tides)  of  0.69  was 
used  (Schwiderski  1980,  Hendershott  1981).  Zero  normal  flow  specified 
boundary  conditions  were  applied  to  all  coastal  and  island  boundaries. 
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The  fully  nonlinear  form  of  the  shallow-water  equation  (Equations  1-3) 
was  used  for  all  grids.  The  nonlinear  bottom  friction  coefficient  Cf  was 
specified  equal  to  0.003  throughout  the  domain.  The  GWCE  parameter 
Xq  was  specified  equal  to  0.001,  which  represents  a  balance  between  the 
primitive  continuity  and  wave  equation  portions  of  the  GWCE  equation. 

Time-steps  used  with  each  grid  are  listed  in  Table  2.  These  time-steps 
were  selected  to  maintain  the  maximum  Courant  number  based  on  wave 
celerity  less  than  unity.  Each  run  was  spun  up  from  a  cold  start  applying  a 
very  smooth  hyperbolic  tangent  ramp  function,  which  acts  over  approxi¬ 
mately  12  days,  to  both  the  boundary  and  direct  forcing  functions.  The 
use  of  this  ramp  function  avoids  exciting  short  period  gravity  modes  and 
vortex  modes  in  the  subinertial  frequency  range  due  to  start-up  transients 
(Reid  and  Whitaker  1981).  Nonetheless,  free  Helmholtz  modes  did  appear 
to  be  excited  in  the  Gulf  of  Mexico  and  the  Caribbean  Sea.  However,  the 
use  of  the  12-day  ramp  minimized  the  amplitude  of  these  modes  and  they 
were  almost  entirely  dissipated  after  about  30  days.  Simulations  were  run 
for  50  days  and  tidal  elevations  were  recorded  at  the  147  defined  elevation 
stations.  The  last  10  days  of  these  tidal  elevation  records  were  analyzed 
using  a  least  square  harmonic  analysis  and  decomposed  into  amplitudes 
and  phases  for  the  M2  tide  and  its  overtides. 

Comparisons  of  M2  elevation  amplitudes  and  phases  are  the  basis  for 
all  grid  intercomparisons.  Global  and  regional  error  measures  established 
by  comparing  the  computed  amplitudes  and  phases  to  a  converged  solu¬ 
tion  (which  for  our  studies  was  obtained  using  grid  R5)  are  presented  in 
Tables  3  and  4.  M2  amplitude  errors  were  computed  as  a  proportional  stan¬ 
dard  deviation  for  each  grid  as  follows: 


pgrid-R5  _ 
11  M2~  amp  ~ 


(11) 


where 


L 


the  total  number  of  elevation  stations  within  a  given 
region 

the  location  of  an  elevation  station 

the  M2  elevation  amplitude  for  a  given  grid  at  station 
coordinates  (xt ,  yx ) 
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(x^.y^)  =  the  M2  elevation  amplitude  for  the  converged  solution 


obtained  using  grid  R5  at  station  coordinates  (xt ,  yx ) 


Eft[ld-amp  may  ^  interpreted  as  the  standard  deviation  in  terms  of  a  frac¬ 
tion  of  a  regional  representative  measure  of  Af2  amplitude  or  roughly  as  an 
average  percentage  of  error  (when  multiplied  by  100)  for  the  M2  ampli¬ 
tude  within  a  specified  region.  M2  phase  errors  were  computed  for  each 
grid  as  an  absolute  average  error  over  a  defined  region  as  follows: 


L 

Egrid-RS  =  1  V 
M2  -phase  ~  £  t-t 
i=l 


(12) 


where 


,yt )  =  the  M2  elevation  phase  for  a  given  grid  at  measurement 
location  (^ ,  y; ) 

(p^5  (j^  ,  3^ )  =  the  Af2  elevation  phase  for  the  converged  solution 

obtained  with  grid  R5  at  measurement  location^ ,  yx ) 


Eft[ld-amp  and  E\ Uld-pRhase  are  Presenled  in  Tables  3  and  4  for  all  grids  con¬ 
sidered  for  the  entire  domain  (i.e.,  at  all  147  elevation  stations  discussed 
in  the  previous  section)  as  well  as  for  various  groups  of  stations  defined  in 
the  previous  section. 


Resolution  Requirements  for  Tidal 
Computations 

In  the  first  numerical  experiment,  the  sensitivity  of  the  solution  to  the 
level  of  refinement  was  examined  using  a  regular  grid  structure  by  com¬ 
paring  simulations  for  grids  R1  through  R5.  Grid  Rl,  shown  in  Figure  3, 
is  very  coarsely  discretized  with  a  resolution  of  approximately  1.6°  x  1.6° 
(this  corresponds  roughly  to  140  km  x  140  km).  The  degree  of  resolution 
for  both  bathymetry  and  the  coastal  boundary  is  identical  to  that  of  the 
grid  itself.  The  only  islands  included  in  grid  Rl  are  Cuba  and  Hispaniola. 
Figure  8  shows  the  A/2  wavelength-to-grid  size  ratio  (  /  Ax )  based  on 

inviscid  linear  long-wave  theory  and  indicates  that  large  portions  of  the 
continental  shelf  are  significantly  underresolved  (regions  where 

Ax<  12.5).  In  fact.  Table  2  indicates  that  the  shallowest  regions  of 

the  domain  are  quite  poorly  resolved,  with  X^/Ax  as  low  as  1.5. 
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Grid  R2  is  derived  from  grid  R1  by  splitting  each  element  into  four  ele¬ 
ments.  This  halves  the  resolution  to  approximately  0.8°  x  0.8°  (roughly 
70  km  x  70  km).  The  coastal  boundary  remains  identical  to  that  of 
grid  R1  although  the  bathymetric  resolution  was  updated  to  be  consistent 
with  the  resolution  in  the  grid.  This  effectively  steepens  the  gradient  of 
the  continental  slope  and  improves  the  representation  of  the  continental 
shelf.  Grids  R3  and  R4  were  subsequently  generated  in  the  same  manner 
as  grid  R2,  each  increasing  the  grid  resolution  by  a  factor  of  2  as  com¬ 
pared  to  the  previous  grid,  while  maintaining  the  same  coastal  boundary 
as  grid  R1  and  updating  the  bathymetry  to  the  resolution  of  the  grid  itself. 
Consequently  grid  R4  is  finely  discretized  over  the  domain  with  a  resolu¬ 
tion  of  approximately  12'  x  12'  (roughly  17  km  x  17  km).  Finally,  grid  R5 
was  generated  by  increasing  the  grid  resolution  by  a  factor  of  2  in  regions 
with  bathymetry  less  than  3,500  m  while  maintaining  the  same  resolution 
as  grid  R4  in  waters  deeper  than  3,500  m.  The  resulting  grid  is  shown  in 
Figure  4.  Figure  9  plots  the  /  Ax  ratio  for  grid  R5.  The  minimum 

Ax  ratio  is  equal  to  approximately  23,  which  should  be  sufficient  to 

resolve  M2  waves  accurately  in  regions  where  topography  does  not  vary 
rapidly.  This  assertion  is  based  on  the  results  of  Fourier  analysis  briefly 
discussed  in  the  section  on  the  governing  equations  and  the  numerical  dis¬ 
cretization.  Figure  9  also  indicates  that  all  regions  in  grid  R5  with  rapidly 
varying  topography,  i.e.,  the  continental  shelf  break  and  slope,  as  well  as 
the  deep  ocean,  are  resolved  with  /  Ax  ratios  ranging  from  approxi¬ 
mately  100  to  over  500. 

The  computed  M2  response  at  the  147  established  elevation  stations 
indicated  that  amplitudes  typically  tend  to  increase  when  comparing 
grids  R1  and  R2  and  may  increase  or  decrease  when  comparing  grids  R2 
and  R3,  grids  R3  and  R4,  and  grids  R4  and  R5.  There  are  a  variety  of 
factors  which  come  into  play  when  considering  the  convergence  of  these 
solutions.  First,  as  the  grid  resolution  improves,  the  M2  wave  is  better  re¬ 
solved  in  the  shallowest  portions  of  the  domain,  at  least  on  a  local  basis, 
and  in  regions  of  rapidly  varying  flow  such  as  in  the  vicinity  of  the  conti¬ 
nental  shelf  break  and  slope.  This  improved  resolution  leads  to  reduced 
truncation  errors  in  these  regions.  Second,  as  grid  resolution  increases, 
the  continental  shelf  and  slope  are  more  accurately  represented.  There¬ 
fore,  the  interaction  between  the  tides  and  topography,  as  well  as  the  im¬ 
portance  of  the  bottom  friction  term,  will  be  altered. 

Error  measures  for  M2  amplitude  and  phase  are  presented  in  Tables  3 
and  4.  In  general,  the  differences  between  solutions  for  subsequent  grids 
become  increasingly  smaller  with  the  finer  grids.  Errors  tend  to  be  largest 
in  the  Gulf  of  Mexico  and  Caribbean  Sea  and  smallest  in  deep  water.  The 
relatively  large  amplitude  and  phase  errors  in  the  Gulf  of  Mexico  are  pre¬ 
dominantly  caused  by  large  errors  at  a  few  stations  along  the  coast  of 
northern  Florida  and  the  Mississippi  delta  (stations  24,  25,  26,  27,  28,  96 
and  98)  in  addition  to  stations  near  the  center  of  the  Gulf  (stations  37,  82, 
and  102).  The  absolute  phase  differences  between  grids  R4  and  R5  associ¬ 
ated  with  these  stations  range  between  8.8  deg  and  17.6  deg.  These  stations 
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lie  in  the  vicinity  of  very  high  gradients  in  phase  which  occur  on  the  Flor¬ 
ida  shelf,  off  the  Mississippi  delta,  and  near  the  center  of  the  Gulf  associ¬ 
ated  with  amphidromes.  The  increased  error  levels  within  the  Caribbean 
are  associated  with  stations  in  the  vicinity  of  Puerto  Rico  and  Hispaniola 
(stations  50  and  54),  which  are  again  located  in  the  vicinity  of  high  gradi¬ 
ents  in  phase  associated  with  an  amphidrome.  Based  on  the  errors  in  Ta¬ 
bles  3  and  4,  the  resolution  provided  in  grid  R5  is  judged  to  allow  M 2 
amplitudes  to  be  predicted  with  less  than  1  deg  error  and  M2  phase  with 
less  than  2  deg  error  (based  on  the  diminishing  error  trends,  since  the  accu¬ 
racy  of  grid  R5  must  be  better  than  that  of  grid  R4)  except  at  a  few  spe¬ 
cific  locations  in  the  Gulf  of  Mexico  and  the  Caribbean  mentioned  above. 
The  M2  constituent  solution  obtained  with  grid  R5  is  considered  to  be  a 
converged  solution  for  most  of  the  WNAT  domain  with  a  coarse  coastal 
boundary.  Finally,  it  is  noted  that  the  errors  in  the  computed  overtides 
were  not  as  good  as  for  the  astronomical  M2  tide.  For  example,  the  Af6 
constituent  differed  by  an  average  of  21  percent  in  amplitude  and  18  deg 
in  phase  between  grids  R4  and  R5  within  the  domain.  Thus  grid  R5  does 
not  lead  to  converged  overtide  solutions. 


Development  of  an  Optimal  Graded  Mesh 

In  the  second  experiment,  optimally  graded  grids  were  developed  with 
the  goal  of  obtaining  M2  solutions  that  were  comparable  to  the  finest  regu¬ 
lar  mesh,  grid  R5.  While  grid  R5  adequately  resolves  the  M2  tide  in  the 
shallowest  portions  of  the  domain,  it  significantly  over-resolves  the  M2 
tide  in  the  deepest  portions  of  the  domain.  In  fact,  Table  2  indicates  that 
there  are  more  than  1 ,322  nodes  per  M2  wavelength  for  the  deepest  waters. 
Over-resolving  a  grid  to  this  extent  is,  of  course,  extremely  inefficient  from  a 
computational  perspective.  Therefore  grids  VI  and  V2  were  developed. 

Grid  VI,  shown  in  Figure  5,  was  generated  by  starting  with  grid  R1 
and  systematically  refining  regions  for  which  X^2  /  Ax  <  25  until  the 

shallowest  regions  within  the  domain  as  well  as  the  coastal  boundary  were 
discretized  to  approximately  the  same  level  of  detail  as  grid  R5.  Thus, 
grid  VI  is  coarsely  resolved  in  deep  water  (1.6°  resolution)  and  very  finely 
resolved  in  the  shallowest  portions  of  the  domain  (5'  resolution).  Regions 
with  steep  gradients  in  bathymetry  such  as  in  the  vicinity  of  the  continen¬ 
tal  shelf  break  and  over  the  continental  slope  were  generally  refined  to 
satisfy  the  minimum  criterion  50  5  X^/  Ax  <  100  .  This  criterion  was 

based  on  one-dimensional  numerical  experiments  and  truncation  error 
analysis,  which  indicated  that  the  X^/ Ax  =  25  is  not  adequate  to  give 

accurate  solutions  over  the  shelf  break  and  slope  (Westerink,  Muccino, 
and  Luettich  1992).  Transition  regions  between  small  and  large  elements 
were  provided  with  additional  resolution  to  create  a  smooth  transition  be¬ 
tween  element  sizes  and  to  control  element  skewness.  Grid  VI  maintains 
the  same  coarse  coastal  boundary  as  grids  R1  through  R5.  Figure  10 
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shows  the  distribution  of  X^/  Ax  contours  for  grid  VI  and  indicates  that 
the  /  Ax  ratio  is  generally  between  25  and  50  on  the  continental  shelf 

and  ranges  between  75  and  100  in  deep  waters.  Over  the  continental  shelf 
break  and  slope,  \m2  1  Ax  generally  varies  between  50  and  100,  although 

in  the  Caribbean  Sea,  the  ratio  is  as  high  as  250. 

The  computed  errors  in  Tables  3  and  4  indicate  that  grid  VI  gives  al¬ 
most  the  same  response  as  grid  R5.  In  fact,  differences  in  A#2  amplitude 
between  grids  R5  and  VI  are  on  the  order  of  2  percent  and  differences  in 
phase  are  on  the  order  of  3  deg.  As  was  the  case  in  the  regular  grid  con¬ 
vergence  studies,  the  Gulf  of  Mexico  and  Caribbean  exhibit  larger  ampli¬ 
tude  and  phase  errors  than  other  portions  of  the  domain.  In  particular,  the 
Gulf  of  Mexico  has  large  phase  errors  associated  with  stations  along  the 
coast  of  northern  Florida  and  the  Mississippi  delta  (stations  24  through  28 
and  96),  in  the  center  of  the  Gulf  (station  37),  in  addition  to  the  southwest¬ 
ern  portion  of  the  Gulf  (stations  82,  99,  102,  127).  Again,  these  errors  cor¬ 
respond  to  the  large  gradients  in  phase  associated  with  amphidromes  on 
the  Florida  shelf  and  in  the  center  of  the  Gulf. 

Finally,  a  second  graded  grid.  Figure  6,  was  developed  to  determine 
whether  insufficient  resolution  over  the  shelf  break,  slope,  and  selected 
deep  regions  was  a  limiting  factor  in  the  agreement  between  grids  R5 
and  VI.  Grid  V2  was  developed  from  grid  VI  by  providing  resolution  of 
at  least  5'  in  the  vicinity  of  the  continental  shelf  break  (typically  between 
100  m  and  500  m)  and  of  at  least  10'  on  the  slope  (typically  between 
500  m  and  4,000  m).  Again,  additional  resolution  was  provided  in  adja¬ 
cent  regions  in  order  to  control  the  ratio  of  adjacent  element  sizes  as  well 
as  element  skewness.  Additional  resolution  was  also  provided  in  the  deep 
regions  of  the  Gulf  and  Caribbean.  Grid  V2  also  maintains  a  coarse 
coastal  boundary.  Figure  1 1  presents  the  Xj#2  /  Ax  ratios  for  grid  V2  and 

indicates  that  the  Ax  ratio  in  the  vicinity  of  the  continental  shelf 

break,  slope,  and  in  deep  waters  in  the  Gulf  and  Caribbean  is  much  greater 
than  for  grid  VI  (i.e.,  X*^/  Ax  is  equal  to  at  least  250  and  often  greater 

than  500),  whereas  other  regions  are  similarly  discretized. 

Tables  3  and  4  compare  the  Af2  amplitude  response  obtained  with 
grids  R5  and  V2  at  the  147  elevation  stations.  The  response  of  the  graded 
grid  V2  is  almost  identical  to  the  response  of  the  graded  grid  VI.  Thus,  the 
resolution  provided  by  grid  VI  over  the  shelf  break,  slope,  and  in  deep 
waters  in  the  Gulf  and  Caribbean  does  not  limit  the  accuracy  of  the  solution. 
In  fact,  it  can  be  concluded  that  grid  V2  unnecessarily  over-resolves  the 
shelf  break  and  slope  regions. 
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Influence  of  the  Resolution  of  the  Coastal 
Boundary  on  the  Afc  Response 


The  third  experiment  examined  the  sensitivity  of  the  solution  to  the  de¬ 
gree  of  resolution  of  the  coastal  boundary  itself.  Grids  R1  through  R5,  as 
well  as  grids  VI  and  V2,  all  represented  the  coastal  boundary  with  the  res¬ 
olution  of  the  coarsest  grid  R1  Therefore,  grid  V3,  shown  in  Figure  7, 
was  generated  from  grid  VI  by  updating  all  coastal  boundaries  to  a  resolu¬ 
tion  of  roughly  5'  (equivalent  to  the  resolution  of  the  most  resolved  re¬ 
gions  in  grid  VI).  The  islands  of  Puerto  Rico,  Andros,  and  Jamaica,  in 
addition  to  50  smaller  islands,  were  also  added.  Tables  3  and  4  indicate 
that  the  elevation  response  changes  substantially  in  some  regions  with 
the  improved  representation  of  the  coastal  boundary.  Particularly  in  near¬ 
coastal  regions  and  on  the  continental  shelf,  the  resolution  of  the  coastal 
boundary  will  strongly  influence  the  computed  response.  Differences  in 
phase  in  the  Gulf  of  Mexico  are  particularly  large.  Differences  in  ampli¬ 
tude  and  phase  are  less  in  remote  and  deep  waters. 
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Comparison  Between 
Computed  and  Observed 
Tidal  Elevations  in  the 
WNAT  Model 


Computations  are  now  presented  using  grid  V3  driven  on  the  open 
ocean  boundary  using  results  from  an  existing  global  model  and  within  the 
interior  domain  with  eight  diurnal  and  semidiurnal  astronomical  tidal  con¬ 
stituents.  These  computations  allow  for  realistic  comparisons  between  the 
coupled  global/WNAT  model  results  and  the  measurement  data  at  the  77 
stations  previously  defined,  since  the  constituents  are  allowed  to  fully  in¬ 
teract  through  the  various  nonlinear  terms  in  the  shallow-water  equations 
(Westerink  et  al.  1989).  The  computations  were  intended  to  be  entirely 
predictive  and  therefore  no  calibration/tuning  procedures  were  performed. 

The  open  ocean  boundary  was  forced  using  tidal  elevations  for  the  K j, 
Ox,  P\,  Q\,  M2,  S2 ,  N2,  and  K2  constituents  from  Schwiderski’s  global 
model  (Schwiderski  1979;  1981,  a-g).  An  effective  tidal  potential  forcing 
within  the  interior  domain  was  applied  for  the  same  eight  constituents  that 
were  forced  on  the  open  ocean  boundary.  Tidal  potential  amplitudes  and 
the  associated  effective  earth  elasticity  factors  for  the  eight  constituents 
are  listed  in  Table  5.  In  particular,  the  diurnal  tides  have  associated  effec¬ 
tive  earth  elasticity  factors  ranging  between  0.695  and  0.782  (Wahr  1981, 
Woodworth  1990)  instead  of  the  widely  used  value  of  0.69  (Schwiderski 
1980,  Hendershott  1981).  Equilibrium  tides  were  directly  computed  and 
therefore  no  nodal  factors  or  astronomical  arguments  were  used  in  either 
the  boundary  or  the  interior  domain  forcing  functions. 

A  constant  bottom  friction  coefficient  equal  to  Cy  =  0.003  was  applied 
throughout  the  domain.  Furthermore,  all  nonlinearities  in  the  shallow-water 
equations  (Equations  1-3)  were  included  in  the  computations.  The  model 
was  spun  up  from  homogeneous  initial  conditions  using  a  12-day  hyperbolic 
tangent  ramp  function.  Due  to  the  presence  of  free  Helmholtz  modes  within 
the  Gulf  of  Mexico  and  Caribbean,  results  were  only  recorded  after  40  days 
when  these  free  modes  had  been  entirely  dissipated.  The  simulations  were 
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run  for  225  days;  during  the  last  185  days  data  were  recorded  every  half 
hour  at  elevation  stations  fcr  which  measurement  data  are  available  (sta¬ 
tions  1-77)  and  every  hour  at  all  nodes  within  the  domain.  A  time-step  of 
37.5  sec  was  used,  giving  a  Courant  number  between  0.01  and  1.1. 

The  computed  time  histories  at  all  nodes  within  the  domain  as  well  as 
at  the  77  measurement  tidal  elevation  stations  were  harmonically  analyzed 
over  the  recorded  185-day  period.  Cotidai  charts  for  the  four  diurnal  constit¬ 
uents  and  the  four  semidiurnal  constituents  are  presented  in  Figures  12 
and  13,  respectively.  These  figures  also  indicate  the  relative  error  between 
the  predicted  and  measured  values  at  the  77  measurement  stations  for 
which  data  are  available.  Figures  12  and  13  indicate  a  marked  similarity 
between  each  of  the  diurnal  and  each  of  the  semidiurnal  constituents.  The 
diurnal  constituent  amplitudes  tend  to  increase  significantly  within  the 
Gulf  of  Mexico.  The  diurnal  constituents  exhibit  an  amphidrome  in  the 
North  Atlantic  off  the  coast  of  Nova  Scotia  as  well  as  dual  amphidromes 
in  the  vicinity  of  the  Bahamas.  The  semidiurnal  constituents  experience 
rapid  increases  in  amplitude  in  the  Gulf  of  Maine,  on  the  Blake  Plateau, 
on  the  Florida  Shelf,  and  on  the  Nicaragua  Rise,  as  well  as  in  the  vicinity 
of  Trinidad.  Semidiurnal  amplitudes  generally  tend  to  be  significantly 
smaller  in  the  Gulf  of  Mexico  and  the  Caribbean  than  in  the  Atlantic.  Full 
or  degenerate  amphidromes  exist  off  of  Puerto  Rico  and  on  the  Florida 
Shelf.  A  full  amphidrome  exists  in  the  center  of  the  Gulf  of  Mexico  for 
the  M2  and  N2  constituents,  while  this  feature  appears  as  a  degenerate 
amphidrome  off  of  Cancun  for  the  S2  and  K2  constituents.  A  full 
amphidrome  exists  in  the  northwest  Gulf  of  Mexico  for  the  K2  constituent. 

The  accuracy  of  the  predicted  tides  was  quantified  by  comparing  pre¬ 
dicted  values  at  the  77  measurement  stations  for  the  eight  astronomical 
constituents  to  the  available  field  data  in  harmonic  form.  The  amplitude 
error  for  constituent  j  for  each  of  the  eight  astronomical  constituents  was 
computed  as  a  proportional  standard  deviation  for  a  defined  region  as: 


tp  V3  -  meas 
j- amp 


-J  [/sj  II 

n  jyj  r- 

T\™as{xx.y J 

2 

(13) 


where 


L  =  total  number  of  elevation  stations  within  a  given  region 
(xl,yl)  =  location  of  an  elevation  station 
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i\y*  (jcx  ,  yx )  =  computed  elevation  amplitude  for  constituent  j  at 
station  coordinates  (jct ,  yt ) 

A 

T\jneas  ^  ^  )  =  measured  elevation  amplitude  for  constituent  j  at 
station  coordinates  (xx ,  yv ) 

The  phase  error  for  each  constituent  j  was  computed  as  an  absolute  aver¬ 
age  error  defined  over  a  region  as: 


p  V3  -meas 
‘-'j  -phase 


i=l 


where 

tp  y3  (xt ,  yt )  =  computed  elevation  phase  for  constituent  j  at 
measurement  location  (jcx  ,  yt ) 

fpjneas  ^  ^  )  =  measured  elevation  phase  for  constituent  j  at 
measurement  location  (jc ,  _yt ) 


(14) 


Amplitude  and  phase  errors  are  presented  in  Tables  6  and  7  for  the  en¬ 
tire  domain  as  well  as  on  a  regional  basis  for  each  of  the  eight  astronomi¬ 
cal  constituents.  The  constituents  are  predicted  with  an  average  amplitude 
error  between  18.2  and  45.3  percent  and  an  average  phase  error  between 
8.3  and  27.5  deg.  Typically,  the  larger  constituents,  the  K j,  O j,  M2,  S2 
and  N2  tides,  have  smaller  relative  amplitude  errors  while  the  smaller  con¬ 
stituents,  the  P],  Q\  and  K2  tides,  have  larger  relative  amplitude  errors. 
There  are  no  predominant  trends  with  respect  to  spatial  variability  in  am¬ 
plitude  errors,  with  the  exception  of  relatively  lower  errors  in  the  Carib¬ 
bean  for  diurnal  constituents  and  relatively  higher  errors  in  the  Caribbean 
and  Gulf  of  Mexico  for  semidiurnal  constituents.  The  distribution  of  am¬ 
plitude  ratio  errors  is  presented  in  Table  8  and  indicates  that  the  diurnal 
constituent  amplitude  predictions  are  evenly  and  closely  clustered  about 
the  measured  values,  while  the  semidiurnal  constituent  amplitude  predic¬ 
tions  are  typically  over-predicted  with  quite  a  few  outliers.  Phase  errors 
for  the  diurnal  tides,  ranging  between  8.6  and  12.4  deg,  are  much  smaller 
than  phase  errors  for  the  semidiurnal  tides,  which  range  between  19.6  and 
27.5  deg.  Phase  errors  tend  to  be  especially  large  for  the  semidiurnal  con¬ 
stituents  within  the  Gulf  of  Mexico  and  Caribbean.  A  consistent  feature 
of  the  results  is  that  stations  which  had  poor  convergence  properties  also 
showed  the  poorest  comparison  to  the  measurement  data.  Again  these 
stations  typically  lie  in  the  vicinity  of  amphidromes.  Finally,  the  relative 
distribution  of  phase  errors  is  presented  in  Table  9,  which  indicates  that 
predicted  diurnal  constituent  phases  generally  lead  and  are  clustered 
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closely  about  the  measured  phase,  while  predicted  semidiurnal  constituent 
errors  are  distributed  over  a  much  wider  range  of  values. 

Resynthesized  time  histories  for  surface  elevation  from  harmonically 
decomposed  computed  constituent  data  (i.e.  resynthesized  ADCIRC-2DDI 
results)  are  compared  to  synthesized  time  histories  of  the  available  mea¬ 
surement  data  at  every  other  station  for  stations  1  through  77  in  Figure  14. 
These  time  histories  were  all  synthesized  using  up  to  eight  astronomical 
constituents  applied  in  the  forcing  functions  or  a  sub-set  of  these  constitu¬ 
ents  when  all  eight  were  not  available  in  the  measurement  data  at  the  sta¬ 
tion.  Procedures  for  resynthesis  are  presented  in  Appendix  A.  Errors  in 
these  time  histories  are  dominated  by  the  constituent  errors  of  the  largest 
constituent  or  constituents.  In  regions  where  semidiurnal  tides  dominate, 
such  as  the  U.S.  east  coast,  the  ADCIRC-2DDI  results  are  almost  always 
over-predicted.  This  relates  directly  to  the  over-prediction  of  the  Af2  con¬ 
stituent.  On  the  other  hand,  in  regions  of  mixed  or  diurnal  tides,  such  as 
in  the  Gulf  of  Mexico,  errors  in  the  time  histories  are  closer  to  the  more 
modest  errors  associated  with  the  dominant  Oj  and  ATj  constituents  in  the 
region. 

On  a  Cray  YMP-6128,  ADCIRC-2DDI  ran  at  0.46  central  processing 
unit  seconds  per  time-step  for  the  WNAT  domain  using  grid  V3. 
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The  overall  accuracy  of  a  numerical  model  depends  on  the  physics  in¬ 
cluded  in  the  model,  on  the  numerical  accuracy  (which  depends  on  the  al¬ 
gorithm  as  well  as  the  grid)  and  on  the  precision  of  the  forcing  functions. 
Two-dimensional  barotropic  tidal  models  generally  give  comparisons  be¬ 
tween  model  predictions  and  measurements  within  10  percent  accuracy 
for  amplitude  and  10  percent  accuracy  for  phase.  Therefore,  it  is  assumed 
that  the  physics  incorporated  in  the  model  used  in  this  report  should  be 
generally  satisfactory. 

The  WNAT  model  grid  convergence  studies  demonstrate  that  conver¬ 
gence  to  solutions  within  1-2  percent  error  in  amplitude  and  2°  to  3°  in 
phase  have  generally  been  achieved,  with  the  exception  of  the  Gulf  of 
Mexico  and  Caribbean,  where  persistent  errors  occur  in  semidiurnal  con¬ 
stituent  phases.  The  latter  problem  appears  to  be  associated  with 
amphidromes,  whose  locations  shift  readily.  The  performance  characteris¬ 
tics  of  the  finest  regular  grid  (6'  x  6'  to  12'  x  12'  in  very  deep  water)  and  the 
graded  VI  and  V2  grids  were  very  similar.  All  three  grids  appear  to  have  dif¬ 
ficulties  in  phase  in  the  Gulf  of  Mexico  and  Caribbean,  with  an  average  error 
of  up  to  10°.  The  results  from  grids  VI,  V2,  and  R5  indicate  that  the  lack  of 
convergence  in  the  Gulf  as  well  as  in  the  Caribbean  is  not  associated  with  in¬ 
sufficient  resolution  over  the  shelf  break,  slope,  in  deep  waters,  or  in  the  im¬ 
mediate  vicinity  of  the  amphidromes.  Rather,  these  problems  appear  to  be 
due  to  factors  that  control  the  structure  and  location  of  the  amphidromes. 
Previous  studies  in  the  Gulf  indicate  that  the  semidiurnal  response  was 
very  sensitive  to  the  grid  resolution  and  bathymetry,  as  well  as  the  bottom 
friction  coefficient  (Westerink  et  al.  1992).  Insufficient  resolution  over 
the  shelf,  particularly  in  the  shallowest  regions,  may  be  responsible  for 
the  relatively  poor  numerical  convergence  in  these  regions.  Finally,  grid 
V3  clearly  demonstrated  the  importance  of  providing  detailed  resolution 
for  the  coastal  boundary.  Phase  response  in  the  Gulf  was  again  the  most 
sensitive. 

Results  of  the  uncalibrated  simulation  in  which  grid  V3  is  forced  with 
Schwiderski’s  (1979,  1 98 1  a-g)  global  model  indicate  that  the  differences 
between  the  WNAT  model  predictions  and  field  data  significantly  exceed  the 
level  of  estimated  numerical  convergence  error  throughout  the  entire  domain. 
Nonetheless  the  largest  errors  again  occur  in  the  Gulf  of  Mexico  as  well  as 
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the  Caribbean.  Many  of  the  largest  semidiurnal  constituent  station  errors 
occur  in  the  vicinity  of  amphidromic  points.  This  suggests  that  the 
amphidromes  are  sensitive  not  only  to  numerical  truncation  errors  and  the 
resolution  provided  for  the  coastal  boundary  but  also  to  the  precision  of 
the  forcing  function.  Indeed,  a  systematic  error  for  all  constituents  exists 
in  amplitude  and  phase  throughout  the  entire  domain  in  deep  as  well  as 
shallow  waters.  In  particular,  all  the  semidiurnal  tides  and  especially  the 
A/2  constituent  are  predominantly  over-predicted.  Thus,  it  is  highly  likely 
that  Schwiderski’s  model  is  not  sufficiently  accurate  on  the  WNAT  model 
open  ocean  boundary.  This  is  not  entirely  surprising,  considering  that 
Schwiderski’s  model  does  not  properly  account  for  dissipation  on  the 
shelf  due  to  a  lack  of  grid  resolution  and  the  linear  character  of  the  model. 
In  fact,  Cartwright  and  Ray  (1991)  computed  root  mean  square  errors  at 
80  deep-ocean  stations  throughout  the  world  and  found  average  amplitude 
errors  of  4.1  cm  for  Af2,  1.86  cm  for  S2,  1.49  cm  for  0\ ,  and  1.76  cm  for 
K\.  How  these  globally  average  errors  exactly  impact  the  constituent  val¬ 
ues  used  on  the  WNAT  model  open  ocean  boundary  is  difficult  to  deter¬ 
mine.  However,  since  the  tides  are  typically  small  in  deep  waters,  errors 
in  Schwiderski’s  results  on  this  boundary  may  be  significant.  Adjusting 
the  constituents  on  the  open  ocean  forcing  boundaries  downwards  from 
Schwiderski’s  values  will  decrease  the  errors  at  the  77  measurement  sta¬ 
tions.  Sensitivity  and  calibration  studies  for  the  forcing  functions  on  the 
open  ocean  boundary  are  currently  being  performed  to  improve  the  match 
between  the  data  and  predictions.  These  findings  will  be  detailed  in  a  sub¬ 
sequent  report. 
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The  WNAT  model  is  a  large  domain  coastal  model  with  a  geometrically 
and  hydrodynamically  simple  open  ocean  boundary.  WNAT  is  ideally  lo¬ 
cated  to  facilitate  tidal  predictions  with  minimum  requirements  for  open 
ocean  boundary  condition  calibration.  The  open  ocean  boundary  lies  al¬ 
most  entirely  in  deep  water  where  astronomical  tides  vary  slowly,  nonlin¬ 
ear  tides  are  small,  and  global  models,  to  which  the  WNAT  model  can  be 
coupled,  are  likely  to  be  most  accurate.  Due  to  the  resulting  large  domain 
and  the  required  resolution  in  shallow  regions  as  well  as  regions  of  rap¬ 
idly  varying  flow,  it  is  highly  desirable  to  vary  nodal  densities  to  fully  re¬ 
solve  flow  features  on  a  local  basis.  The  FE  shallow-water  equation 
model,  ADCIRC-2DDI,  easily  accommodates  the  required  highly  flexible 
meshes. 

Resolution  requirements  have  been  systematically  studied  using  the  en¬ 
tire  WNAT  model  domain  and  a  sequence  of  regular  grids  (with  grids  rang¬ 
ing  from  a  very  coarse  1.6°  x  1.6°  mesh  to  a  very  fine  6'  x  6'  to  12'  x  12' 
mesh)  as  well  as  a  sequence  of  non-uniformly  graded  grids  with  resolution 
varying  between  1.6°  and  5'  within  each  mesh.  Elevation  response  compu¬ 
tations  with  the  sequence  of  regular  grids  indicate  that  the  semidiurnal  A/2 
tide  typically  converges  to  better  than  1  percent  accuracy  in  elevation  am¬ 
plitude  and  2  deg  in  phase  using  the  6'  x  6'  to  12'  x  12'  mesh  with  the  ex¬ 
ception  of  specific  locations  in  the  Gulf  of  Mexico  and  Caribbean  in 
which  shifting  amphidromes  lead  to  consistently  larger  phase  errors.  The 
overtide  constituents  did  not  satisfactorily  converge,  even  with  the  finest 
6'  x  6'  to  12'  x  12'  mesh.  The  elevation  response  computed  using  the  non- 
uniformly  graded  meshes  achieves  about  the  same  level  of  performance  as 
the  finest  nearly  uniform  grid.  Overall,  these  studies  indicate  that  grids 
can  be  developed  based  on  a  X  /  Ax  =  25  criterion  with  some  additional 
resolution  over  the  shelf  break  and  slope.  This  is  in  sharp  contrast  with 
one-dimensional  studies,  which  indicated  the  need  for  a  much  higher  level 
of  resolution  over  the  continental  shelf  break  and  slope  (Westerink, 
Muccino,  and  Luettich  1992b).  This  indicates  the  importance  of  perform¬ 
ing  grid  convergence  studies  in  two  dimensions  as  opposed  to  one  dimen¬ 
sion.  However,  it  appears  that  the  shallow  regions  in  the  Gulf  of  Mexico 
and  Caribbean  need  to  be  more  highly  resolved  in  order  to  converge  to  the 
same  level  of  precision  as  in  other  parts  of  the  domain.  Finally,  it  was 
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noted  that  the  representation  of  the  coastal  boundary  can  also  substantially 
influence  tidal  response. 

Predictive  computations  were  performed  with  the  optimal  graded  mesh 
driven  by  the  eight  constituents  from  Schwiderski’s  global  ocean  model 
results  and  by  tidal  potential  functions.  Comparisons  of  measurement 
data  at  77  stations  for  the  K\,  Oj,  Pj,  Qj,  Af2,  S2,  N2  and  K2  constituents 
indicate  that  these  constituents  are  predicted  with  average  errors  in  ampli¬ 
tude  between  18.2  percent  and  45.3  percent  and  average  errors  in  phase  be¬ 
tween  8.3°  and  27.5°.  Phase  errors  in  the  Gulf  of  Mexico  and  Caribbean 
for  semidiurnal  constituents  are  quite  sensitive  to  the  shifting  amphidromes. 
Even  in  deep  water,  there  is  no  close  match  between  predictions  and  mea¬ 
surements.  It  is  felt  that  Schwiderski’s  model  results  contribute  substan¬ 
tially  to  these  error  levels.  The  sensitivity  of  response  to  perturbations  in 
the  boundary  forcing  function  and  driving  the  model  boundary  with  more 
accurate  global  models  are  currently  being  investigated. 
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Figure  1.  The  Western  North  Atlantic  Tidal  (WNAT)  model  computational  domain, 

including  the  locations  of  tidal  elevation  stations  with  available  field  data  and 
tidal  elevation  stations  without  field  data 


Figure  2.  Topography  (in  meters)  within  the  WNAT  computational  domain 


Figure  3.  Grid  Rl:  Regular  finite  element  1.6°  x  1.6°  discretization  of  the  domain  with  a 
coarse  representation  of  the  coastal  boundary 


Figure  4.  Grid  R5:  Nearly  regular  finite  element  6'  x  6'  to  12'  x  12'  discretization  of  the 
domain  with  a  coarse  representation  of  the  coastal  boundary 


Figure  5.  Grid  VI :  Graded  finite  element  discretization  of  the  domain  based  on  a 

/Ax  2: 25  criterion  with  a  coarse  representation  of  the  coastal  boundary 


Figure  6.  Grid  V2:  Graded  finite  element  discretization  of  the  domain  based  on  a 

/Ax  >25  criterion  with  additional  resolution  in  the  vicinity  of  the  continental 

shelf  break  and  with  coarse  representation  of  the  coastal  boundary 


Figure  7.  Grid  V3:  Graded  finite  element  discretization  of  the  domain  based  on  a 

/Ax>  25  criterion  with  additional  resolution  in  the  vicinity  of  the  continental 

shelf  break  and  with  fine  representation  of  the  coastal  boundary 


Figure  9.  Distribution  of  X^  /Ax>  25  contours  for  grid  R5 


Figure  12.  Computed  contours  for  elevation  amplitudes  (in  meters)  and  phases  (in  degrees 
relative  to  GMT)  for  the  four  diurnal  astronomical  tidal  constituents  used  in  the 
boundary  and  interior  domain  forcing  functions  (Sheet  1  of  8) 
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Figure  12.  (Sheet  7  of  8) 


Figure  12.  (Sheet  8  of  8) 


Figure  13. 


Computed  contours  for  elevation  amplitudes  (in  meters)  and  phases  (in  degrees 
relative  to  GMT)  for  the  four  semidiurnal  astronomical  tidal  constituents  used  in 
the  boundary  and  interior  domain  forcing  functions  (Sheet  1  of  8) 
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Figure  1 3.  (Sheet  3  of  8) 
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Figure  14.  Comparisons  of  time  histories  of  resynthesized  surface  elevations  from  harmon¬ 
ically  decomposed  computed  constituent  data  and  of  synthesized  surface 
elevations  from  harmonically  decomposed  measured  constituent  data  (using 
the  K 01t  P-,,  01 ,  M2,  N2,  S2,  and  Kz  constituents)  at  representative  tidal 
stations  within  the  Western  North  Atlantic  tidal  model  domain  (Sheet  1  of  39) 
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Table  1 

Location  of  Tidal  Elevation  Stations  with  Measured  Data 


Station 

No. 


Station  Name 

Position 

(longitude/ 

latitude) 

Location 

Yarmouth,  Nova  Scotia 

66°7'  -  43°50' 

Coast 

St.  John,  New  Brunswick 

66°3'  -  45°16' 

Coast 

Bar  Harbor,  ME 

68°12'  -  44°24' 

Coast 

Rockland,  ME 

69°6'  -  44°6' 

Coast 

Portsmouth,  Navy  Yard,  NH 

70°44'  -  43°5' 

Coast 

Woods  Hole,  Ocean  Institute,  MA 

70°40'  -  41  °30' 

Coast 

Nantucket,  Nantucket  Island,  MA 

70°6'  -  41°17' 

Coast 

Montauk,  NY 

71°58'  -  41  °3' 

Coast 

Sandy  Hook,  NJ 

74°1'  -  40°28' 

Coast 

Atlantic  City,  NJ 

74°25'  -  39°21' 

Coast 

Cape  May  Ferry  Terminal,  NJ 

74°57'  -  38°58' 

Coast 

Indian  River  Inlet,  DE 

75°4'  -  38°36' 

Coast 

Ocean  City  Fishing  Pier,  MD 

75“5'  -  38°19' 

Coast 

Chesapeake  Bay  Bridge  Tunnel,  VA 

76°7'  -38°58' 

Coast 

Duck  Pier,  NC 

75-45' -36-11' 

Coast 

Southport,  NC 

78°1'  -33°55' 

Coast 

Charleston,  SC 

79-55'  -  32°47' 

Coast 

Savannah  River  Entrance,  GA 

80-54'  -  32°2' 

Coast 

Mayport,  FL 

81-26' -30°24' 

Coast 

Miami  Harbour  Entrance,  FL 

80°8'  -25°46' 

Coast 

Key  West,  FL 

81-48'  -  24-33' 

Coast 

Naples,  FL 

81-48'  -  26°8' 

Coast 

Cedar  Key,  FL 

83°2'  -  29°8' 

Coast 

St.  Marks  Light,  FL 

84°11'-  30°4' 

Coast 

Alligator  Bayou,  FL 

85°45'  -  30-10' 

Coast 

Bay  St.  Louis,  MS 

89°20' -30°18' 

Coast 

Cat  Island,  MS 

89-10'-  30-14' 

Coast 

Southwest  Pass,  LA 

89°25'  -  28°56' 

Coast 

Source1  of  Tidal 
Data 


IHO,  USGS,  NOAA 


IHO,  NOAA 


IHO,  USGS,  NOAA 


IHO,  NOAA 


NOAA 


IHO,  NOAA 


IHO,  USGS,  NOAA 


IHO.  USGS,  NOAA 


IHO.  USGS,  NOAA 


IHO,  NOAA 


USGS.  NOAA 


NOAA 


NOAA 


IHO,  NOAA 


IHO,  NOAA 


IHO,  NOAA 


IHO,  NOAA 


IHO,  NOAA 


IHO.  GOM,  NOAA 


IHO,  GOM 


IHO,  GOM 


IHO,  GOM,  NOAA 


IHO,  GOM 


IHO.  GOM 


IHO,  GOM 


IHO,  GOM 


(Continued) 


1  Sources  of  Tidal  Data.  IHO — Internalional  Hydrographic  Organization  Tidal  Constituent  Bank  (1991); 
USGS:  U.S.  Geological  Survey  (1984);  NOAA:  National  Ocean  Survey,  National  Oceanographic  and 
Atmospheric  Administration;  GOM:  Reid  and  Whitaker  (1981). 


Table  1  (Continued) 


Station 

No. 

Station  Name 

Position 

(longitude/ 

latitude) 

Location 

Source  of  Tidal 

Data 

29 

Point  au  Far,  LA 

91=21'-  29°20' 

Coast 

IHO,  GOM 

30 

Galveston,  TX 

94=47'  -  29°19' 

Coast 

IHO,  GOM,  NOAA 

31 

Port  Aransas,  TX 

97=3'  -  27°49' 

Coast 

GOM 

32 

South  Padre  Island,  TX 

97=9'  -  26=4' 

Coast 

IHO,  GOM 

33 

Ciudad  Madero,  Mexico 

97°51'  -  22°13' 

Coast 

GOM 

34 

Coatzacoalcos,  Mexico 

94°25'  -  18=9' 

Coast 

IHO,  GOM 

35 

Campeche,  Mexico 

90°32'  -  19=50' 

Coast 

IHO,  GOM 

36 

Progreso,  Mexico 

89°40'  -  21=17' 

Coast 

IHO,  GOM 

37 

IAPSO  #30-1 .2.6,  GOM  Pelagic 

89=39'  -  24=46' 

Deep  ocean 

IHO, GOM 

38 

IAPSO  #30-1 .2.13,  Western  Florida  Outer 
Shelf 

84=15'  -  26°42' 

Shelf 

IHO,  GOM 

39 

Puerto  Cortes,  Honduras 

87=57'  -  15°50' 

Coast 

IHO 

40 

Puerto  Cabezas,  Nicaragua 

83=22'  -14=1' 

Coast 

IHO 

1  41 

Puerto  Limon,  Costa  Rica 

83°2'  -9°60' 

Coast 

IHO 

1  42 

Cristobal  (Canal  Zone),  Panama 

79°55'  -  9=21' 

Coast 

IHO 

1  43 

Cartagena,  Colombia 

75°32'  -  10°23' 

Coast 

IHO 

I  44 

Willemstad,  Curacao,  Netherlands  Antilles 

68=56'-  12=6' 

Coast 

IHO 

La  Guaira,  Venezula 

66=56'  -  10°37' 

Coast 

IHO 

46 

Cumana,  Venezuela 

64°1 0'  -  10=27' 

Coast 

IHO 

47 

Port  of  Spain,  Trinidad 

61=31'  -  10=39' 

Coast 

IHO 

48 

Castries,  St.  Lucia,  B.W.I. 

61=0'  -  14°1' 

Coast 

IHO 

49 

Fort-de- France,  Martinique 

61=3'  -  14°35' 

Coast 

IHO 

50 

St.  Thomas,  Virgin  Islands 

64°56'  -  1 8°20' 

Coast 

IHO 

51 

Magueyes  Island,  Puerto  Rico 

67°3'  -  17=58' 

Coast 

IHO 

52 

Ciudad  Trujillo,  Dominican  Republic 

69°53'  -  18=28' 

Coast 

IHO 

53 

Puerto  Plata,  Dominican  Republic 

70=41'  -  19=45' 

Coast 

IHO 

54 

Port-au-Prince,  Haiti 

72=21'-  1 8°33' 

Coast 

IHO 

55 

Guantanamo  Bay,  Cuba 

75=9'  -  19=54' 

Coast 

IHO 

56 

Gibara,  Cuba 

76=7'  -  21 °6' 

Coast 

IHO 

57 

Casilda.Cuba 

79=59'  -  21°45' 

Coast 

IHO 

(Continued) 


Table  1  (Concluded) 


Station 

No.  Station  Name 

58  Havana,  Cuba 

59  Nassau,  Bahamas 

60  Eleuthera,  Bahamas 

61  Ireland  Island,  Bermuda 

62  St.  David’s  Island,  Bermuda 

63  IAPSO  #30-1 .3.2 

IAPSO  #30-1 .2.10,  South  Bermuda  Rise 

64  Region 

65  IAPSO  #30-1 .2.13 

IAPSO  #30-1 .2.14,  Rim  of  Biake 

66  Escarpment 

IAPSO  #30-1 .2.4,  Hatteras  Plain  - 

67  Bermuda  Rise 

68  IAPSO  #30-1 .2.5,  Bermuda  Rise  Region 

IAPSO  #30-1 .2.15,  Rim  of  Blake 

69  Escarpment 

70  IAPSO  #30-1 .2.1 1 ,  East  of  Blake  Plateau 

71  IAPSO  #30-1 .2.18,  Near  Bermuda 

72  IAPSO  #30-1 .2.3 

73  IAPSO  #30-1 .2.16,  East  of  Cape  Hatteras 

74  IAPSO  #30-1 .2.19,  East  of  New  York  Bight 

75  IAPSO  #30-1 .2.17,  Outer  New  York  Bight 

76  IAPSO  #30-1 .2.1 ,  Off  Long  Island,  NY 

IAPSO  #30-1 .2.32 


Position 

(longitude/  Source  of  Tidal 


latitude) 

Location 

82=20'  -  23°9' 

Coast 

77°21'  -  25°5' 

Coast 

76°9'  -  24°56' 

64°50'  -  32°19' 

64°39'  -  32°22' 

E 

64°53'  -  16°32' 

Toe  of  slope 

69=20'  -  26°28' 

Deep  ocean 

84=15'-  26=42' 

Shelf 

Data 


76°47'  -  28=1' 

69°45'  -  28°8' 
67°32'  -  28°14' 


Slope 

IHO 

Deep  ocean 

IHO 

Deep  ocean 

IHO 

76°48'  -  28°27'  Slope  I  HO 


76=25'  -  30=26' 

Toe  of  slope 

IHO 

64=26'  -  32=1' 

Deep  ocean 

IHO 

75=37'-  32=41' 

Toe  of  slope 

IHO 

73=5'  -  37=22' 

Slope 

IHO 

IHO 

72=10'-  39=13' 

Slope 

IHO 

68=38'  -  40=7' 

Shelf  break 

IHO 

77 


70=54'  -  40=1 8' 


Shelf 


IHO 


Table  2 

Properties  of  Grids  R1  through  R5  and  VI  through  V3 


Approximate  Grid 
Size 


Grid 

Nodes 

Boundary 

R1 

434 

Coarse 

R2 

1611 

Coarse 

R3 

6191 

Coarse 

R4 

24255 

Coarse 

R5 

64065 

Coarse 

VI 

10147 

Coarse 

V2 

27816 

Coarse 

V3 

19858 

Fine 

Resc  ition 


Graded 


Graded 


Graded 


5'— >1.6 


5'->1 .6' 


Time- 

step 

min 


Table  3 

M2  Elevation  Amplitude  Errors,  ,  for  the  Grids  Applied  in  the 

Convergence  Study  as  Compared  to  the  Converged  Solution  Obtained 
with  Grid  R5  for  Various  Regions  with  the  Domain 


Entire 

Domain 

Stations 

1-147 


Ml 


Stations  with 
Measurement  Data 


U.S. 

Eas 
Coa 
Stations 
1-20 


Gulf  of 
Mexico 
Stationa 
21-38 


Carib¬ 

bean 

Sea 

Stations 

39-58 


Remote 

Stations 

59-77 


Additional  Stations 
without  Measurement  Data 


Vicinity  Vicinity 

Conti-  of  the  of  the 

nental  Shelf  Toe  of  Deep 

Shelf  Break  Slope  Ocean 

Stations  Stations  Stations  Stations 

78-87  88-109  110-130  131-147 


Table  4 

M2  Elevation  Phase  Errors,  ^u'-phase  » *or  *he  Grids  Applied  in  the 

Convergence  Study  as  Compared  to  the  Converged  Solution  Obtained 
with  Grid  R5  for  Various  Regions  within  the  Domain 


EM^pfhas»  •  degree 


Stations  with 
Measurement  Data 


Additional  Stations 
without  Measurement  Data 


Domain  Coast 


Gulf  of 
Mexico 
Stations 
21-38 

— 

Carib¬ 

bean 

Sea 

Stations 

39-58 

Remote 

Stations 

59-77 

30.1 

29.1 

8.4 

12.9 

21.1 

2.3 

6.8 

5.6 

0.95 

6.3 

3.2 

0.31 

10.4 

2.6 

0.33 

9.6 

1.8 

0.53 

41.5 

14.2 

2.2 

Vicinity 
of  the 
Toe  of 
Slope 


Deep 

Ocean 


Table  5 

Tidal  Potential  Constants  for  the  Principal  Tidal  Constituents 
(from  Reid  (1990))  and  the  Associated  Effective  Earth  Elasticity  Factor 
(from  Wahr  (1981)) 


K2  luni-solar 


Species,  j 

n 

Constituent 

T/« ,  hr 

C/n  ,  m 

0 ijn 

1 

i 

Kf  luni-solar 

23.934470 

0.141565 

0.736 

2 

Of  principal  lunar 

25.819342 

0.100514 

0.695 

3 

Pf  principal  solar 

24.065890 

0.046843 

0.706 

4 

01  elliptical  lunar 

26.868357 

0.019256 

0.695 

2 

1 

Mz  principal  lunar 

12.420601 

0.242334 

0.693 

2 

S2  principal  solar 

12.000000 

0.112841 

0.693 

3 

N2  elliptical  lunar 

12.658348 

0.046398 

0.693 

11.967235 


0.030704 


0.693 


Table  6 

Amplitude  Errors,  ™as ,  for  Constituent  /  Computed  Using  Grid  V3  as 
Compared  to  Field  Measured  Values  for  Various  Regions  within  the  Domain 


_  V3 -meas 
j-amp 

Entire  Domain 
Stations  1-77 

U.S.  East 

Coast 

Stations  1-20 

Gulf  of  Mexico 
Stations  21-38 

Caribbean  Sea 
Stations  39-58 

Remote 

Stations  59-77 

0.182 

0.220 

0.199 

0.077 

0.131 

0.205 

0.212 

0.218 

0.109 

0.199 

0.345 

0.289 

0.453 

0.167 

0.235 

0.324 

0.336 

0.356 

0.131 

0.251 

0.270 

0.266 

0.278 

0.652 

0.204 

0.290 

0.288 

0.407 

0.529 

0.122 

Table  7 

Phase  Errors,  ,  for  Constituent  /Computed  Using  Grid  V3  as 

Compared  to  Field  Measured  Values  for  Various  Regions  within  the  Domain 


U.S.  East 
Coast 

Stations  1-20 


V3-meas 

£/-p/ies« 


Gulf  of  Mexico  Carlbbesn  Sea  Remote 
Stations  21-38  Stations  39-58  Stations  59-77 


8 


23.6 


35.8 


12.7 


Table  8  a  a 

Distribution  of  Amplitude  Ratio  Error,  Rg  =  rjyV3(xl ,  yx)  /tij"eas  (xt ,  yt) ,  for  the  Measurement  Stations 


10.7  7.1  28.6 


Table  9 

Distribution  of  Phase  Error,  Pd  =  (p,w(x  ,  yt)  -  <p,mMS(x  ,  yt) ,  for  the  Measurement  Stations 


Appendix  A 

Resynthesis  of  Time  Histories 
Using  the  Cotidal  Charts 


In  order  to  resynthesize  a  time  history  from  the  cotidal  charts  presented 
in  Figures  12  and  13  in  the  main  text  and  the  more  detailed  cotidal  charts 
for  the  U.S.  east  coast  in  Figures  B1  and  B2  and  the  Gulf  of  Mexico  in  Fig¬ 
ures  Cl  and  C2,  the  following  steps  are  applied. 

Step  1 


Interpolate  amplitudes  and 
phases  for  the  eight  astronomical 
constituents  presented  in  this 
report  using  the  most  detailed 
set  of  cotidal  charts  available. 

For  example,  in  order  to  obtain 
values  for  Duck  Pier,  NC  (lo¬ 
cated  75°7'W  and  36°11'N),  the 
set  of  cotidal  charts  presented 
in  Figures  B1  and  B2  are  inter¬ 
polated.  The  adjacent  ampli¬ 
tudes  A  i  and  phases  T,  can  be 
readily  obtained: 

Step  2 

Surface  elevations  as  a  function  of  time  can  now  be  computed  using 
the  following  relationship: 

8 

C  (*>  y.  t)  =  X  Ai  (*>  y)  fi  (fo)  . . , , 

f  (Al) 

cos  [to,-  (t  -  r0)  +  Vj  (t0)  -  '¥{.rad,ans)  (Xj ,)] 


Constituent 

Ag ,  rn 

M'/.deg 

0.092 

179 

Ot 

0.066 

183 

0.038 

173 

0.012 

179 

n2 

0.14 

354 

m2 

0.63 

14 

s2 

0.11 

57 

k2 

0.037 

48 

Appendix  A  Resynthesis  of  Time  Histories 


Al 


where  £(*,  y,  t)  equals  the  free  sur¬ 
face  elevation  at  the  site  of  interest 
as  a  function  of  time,  A(( x,  y)  equals 
the  amplitude  for  the  ilh  constituent, 

'Fj.(jc,  y)  equals  the  phase  for  the  ith 
constituent  (in  radians  per  second), 

/((r0)  equals  the  nodal  factor  rela¬ 
tive  to  reference  time  /0  for  the  ith 
constituent,  and  v^tg)  equals  the 
equilibrium  argument  relative  to  ref¬ 
erence  time  t0  for  the  ith  constitu¬ 
ent.  For  the  eight  constituents  of 
interest,  0);  has  the  following  values: 

Guidelines  as  to  how  to  compute  ffF  and  v?F  are  given  by  Schureman 
(1941)  and  Foreman  (1977).  If  equilibrium  tides  are  computed,  fFF  and 
should  be  specified  as  1.0  and  0.0  for  all  constituents.  A  resynthe¬ 
sized  signal  for  the  data  is  presented  between  50  and  60  days  in  Figure  Al. 


Constituent 

cof ,  radians/second) 

*1 

0.000072921158358 

o 1 

0.000067597744151 

Pi 

0.000072522945975 

Oi 

0.000064958541129 

n2 

0.000137879699487 

m2 

0.000140518902509 

s2 

0.000145444104333 

K 2 

0.000145842317201 

Step  2  (Alternate) 

An  alternate  procedure  to  computing  a  time  series  according  to  Equa¬ 
tion  Al  is  to  use  the  constituent  tide  record  generation  program  of  the  Au¬ 
tomated  Coastal  Engineering  System  (ACES)  (Leenknecht,  Szuwalski, 
and  Sherlock  1990).  This  program  requires  constituent  data  for  the  partic¬ 
ular  desired  location  to  predict  a  tide  elevation  at  that  locale  for  the  period 
of  time  specified  by  the  user.  The  tidal  constituents  contained  in  this  re¬ 
port  can  be  used  as  input  to  the  ACES  program  with  one  exception.  The 
Greenwich  epochs  shown  in  Equation  Al  must  be  converted  to  local  ep¬ 
ochs  K.  This  transformation  is  according  to  the  following  relationship: 


Greenwich  epoch 'Vj  =  K  +  pL  (A2) 

where  k  is  the  local  epoch  needed  for  input  to  the  ACES  program  to  re¬ 
construct  a  local  tide,  p  is  the  coefficient  of  the  constituent  (i.e.,  2  for  Af2, 

1  for  Kj ,  etc.),  and  L  is  the  longitude  of  the  location  for  which  a  tidal  pre¬ 
diction  is  desired.  By  making  the  above  adjustment,  this  report  can  be 
used  in  conjunction  with  the  ACES  program  to  yield  tidal  predictions  at 
any  location  along  the  east  coast  of  the  United  States  and  Gulf  of  Mexico. 


A2 


Appendix  A  Resynthesis  o<  Time  Histories 


Figure  A1.  Sample  time  history  of  resynthesized  surface  elevations  from  harmonically 

decomposed  computed  constituent  data  (using  the  AC, ,  O, ,  P, ,  O, ,  M2,  N2,  S2, 
and  K2  constituents)  at  an  example  station  within  the  Western  North  Atlantic 
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Appendix  B 

Cotidal  Charts 

for  the  U.S.  East  Coast 


Appendix  B  Cotidal  Charts  for  the  U.S.  East  Coast 


Figure  B1.  Computed  contours  for  elevation  amplitudes  (in  meters)  and  phases  (in  degrees 
relative  to  GMT)  for  the  ^ ,  01 ,  P1 ,  and  Q1  diurnal  astronomical  tidal  constituents 
for  the  U.S.  east  coast  (Sheet  1  of  8) 


B2 
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Figure  B1.  (Sheet  2  of  8) 
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Figure  B1 .  (Sheet  4  of  8) 
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Figure  B1 .  (Sheet  5  of  8) 
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Appendix  B  Cotidal  Charts  for  the  U  S.  East  Coast 


Figure  B1 .  (Sheet  6  of  8) 
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Figure  B1 .  (Sheet  8  of  8) 


Appendix  B  Cotidal  Chads  for  the  U.S.  East  Coast 


Figure  B2.  Computed  contours  for  elevation  amplitudes  (in  meters)  and  phases  (in  degrees 
relative  to  GMT)  for  the  A/1 ,  S1 ,  and  /C,  semidiurnal  astronomical  tidal 
constituents  for  the  U.S.  east  coast  (Sheet  1  of  8) 


BIO 
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Figure  B2.  (Sheet  2  of  8) 
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Figure  B2.  (Sheet  6  of  8) 
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Figure  B2.  (Sheet  8  of  8) 
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Cotidal  Charts 
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Appendix  C  Cotidal  Charts  for  the  Gulf  of  Mexico 


Cl 


Figure  Cl.  Computed  contours  for  elevation  amplitudes  (in  meters)  and  phases  (in  degrees 
relative  to  GMT)  for  the  Kv  01t  P^,  and  Q1  diurnal  astronomical  tidal  constituents 
for  the  Gulf  of  Mexico  (Sheet  1  of  8) 
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Figure  Cl.  (Sheet  2  of  8) 
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Figure  Cl.  (Sheet  3  of  8) 
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Figure  Cl .  (Sheet  6  of  8) 
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Figure  C2.  Computed  contours  for  elevation  amplitudes  (in  meters)  and  phases  (in  degrees 
relative  to  GMT)  for  the  Mv  N^,  St,  and  semidiurnal  astronomical  tidal 
constituents  for  the  Gulf  of  Mexico  (Sheet  1  of  8) 
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